Majorana zero modes in dislocations of Sr 2 Ru0 4 
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We study the topologically protected Majorana zero modes induced by lattice dislocations in chiral topologi- 
cal superconductors. Dislocations provide a new way to realize Majorana zero modes at zero magnetic field. In 
particular, we study several different types of dislocations in the candidate material Sr 2 Ru04. We also discuss 
the properties of linked dislocation lines and linked dislocation and flux lines. Various experimental conse- 
quences are predicted which provide a new approach to determine whether nature of the superconducting phase 
of Sr 2 Ru0 4 . 



Introduction: Chiral topological superconductors (TSC) in 
2D have garnered much attention over the past decade after the 
prediction that vortices in these superconducting phases har- 
bor bound Majorana zero-modes| 1 ]. This property is an im- 
portant feature for topological quantum computing architec- 
tures which are based on the use of such non-Abelian anyon 
qubits(21|3l. It is generally agreed that Sr2Ru04 is a quasi- 
2D, p-wave superconductor with broken time-reversal sym- 
metry, although the precise nature of the order parameter is 
still controversial (for a review, see Ref. (4J). So, while there 
is no concrete evidence that Sr 2 Ru04 is a topological super- 
conductor, it is one of the best candidate chiral topological 
superconductor materials. 

The low-energy electronic structure of the normal metal 
state of Sr 2 Ru04 is controlled by the ti g multiplet of d- 
orbitals d xz , d yz , and d xy . These three orbitals give rise 
to three Fermi-surfaces which are expected to become fully 
gapped below the superconducting transition temperature at 
T c ~ 1.5 K. The conventional wisdom indicates that the 
quasi-2D d xy band dominates the pairing instability and de- 
velops a nodeless chiral p x + ip y order parameterQ. If 
such an order parameter were generated, then the recently 
measured half-quantum vortices would indeed bind Majorana 
fermion states! 6 , 7 ] which could be manipulated as qubits. On 
the other hand, a conflict between the theoretical prediction of 
chiral surface states in the p x + ip y state, and the clear lack 
of surface currents measured in Ref. [8], motivated Raghu, 
Kapitulnik, and Kivelson (RKK) to make a different, com- 
pelling proposal for the nature of the order parameter which 
is dominated instead by the quasi- ID d xz and d yz orbitals (9). 
The RKK order parameter is again, odd-parity/p-wave which 
breaks time-reversal, but it is not a chiral topological super- 
conductor and should not exhibit chiral edge states. Thus, it is 
consistent with the measurements of Ref. 1 8 ]. The true nature 
of the order parameter is still an open question, and in light of 
this one seeks experiments that can distinguish between these 
two theories. 

In this Letter we propose new experiments that can dis- 
tinguish the two pairing schemes and also provide a new 
approach to Majorana zero modes. Our key observation is 
that the RKK order parameter, while trivial in the sense of 
2D chiral superconductor, in fact is still nontrivial as a weak 



topological superconductor. Weak topological superconduc- 
tors (and insulators) (lOHIS are topological states protected 
by translation symmetry. They are distinguished by a topo- 
logical invariant defined in a lower-dimensional sub-manifold 
of the Brillouin zone (BZ) (recall that the so-called strong in- 
variant depends on the electronic structure in the entire BZ). 
For example, the 3D weak topological insulators are classi- 
fied by three Z2 invariants, which characterize whether the 
gapped Bloch Hamiltonian restricted to the three 2D planes 
of k x = 7r, k y = 7r, k z = 7r in the BZ are trivial or quantum 
spin Hall insulators in that plane. Heuristically, a non-trivial 
weak invariant indicates a state which is made from stacking 
up topological states of lower dimensions. Another example 
is a 2D superconductor in class D (no symmetry) which has 
two weak Z2 indices. These exist since there is a strong Z2 
topological superconductor in this class in 1 -dimension -the 
Majorana chain/p-wave wire fT3llT5ll . We show that the RKK 
pairing of d xz and d yz bands corresponds to this type of weak 
topological superconductor. As a consequence, we show that 
naturally occurring or fabricated crystal defects exhibit a num- 
ber of remarkable properties that are not shared by the conven- 
tional p x + ip y state. These properties are not only useful for 
a characterization of the superconducting state but may have 
applications in quantum computing since, as we will show, 
these defects can trap Majorana zero modes. 

3D weak topological superconductors: The non- 
interacting topological insulators and superconductors in 
generic dimensions have been classified (13] [15J [16) . Here we 
are interested in superconductors of class D (without time- 
reversal or spin-rotation symmetry), which have a strong Z2 
topological classification in ID, Z in 2D, and trivial in 3D. For 
class D SCs in 3D with translational symmetries we can de- 
fine weak topological invariants as well — the Chern numbers 
(which are defined in 2D) can be defined along constant k x , k y 
or k z planes in the BZ. In a gapped state the Chern number 
cannot change, so the Chern number in different k z = const 
planes is the same integer n z . Similarly n x and n y can be de- 
fined for the other two planes, as is illustrated in Fig. [IJa). The 
integer- valued vector n = (n x , n yi n z ) are the primary weak 
indices of the 3D TSC. A system with indices n is topologi- 
cally equivalent to a set of decoupled 2D layers of topologi- 
cal chiral superconductors with non-vanishing Chern number, 



FIG. 1. (a) The Brillouin zone of a 3D SC. The primary weak topo- 
logical invariants n = (n x ,n y ,n z ) are defined as the Chern numbers 
in the three independent planes (colored green, blue and red), and 
the secondary weak topological invariants v — (is x ,is y ,is z ) are de- 
fined as the ID Z 2 indices along the three perpendicular lines colored 
in red. (b) Illustration that a topological superconductor with topo- 
logical invariants n = (0,0, 1), v — (1, 1,0), such as Sr2Ru04 
(see text) is equivalent to decoupled layers of chiral superconduc- 
tors (with the red line and arrow labeling the chiral edge states) and 
ID wires along x and y directions with Majorana zero modes at end 
points. 



stacked along the n direction. For any surface plane which is 
not perpendicular to n, there will be chiral surface states. 

Similarly, the ID Z 2 invariants fT4l can be calculated along 
time-reversal invariant lines in the 3D BZ[ 17]. The three sec- 
ondary weak topological invariants are defined as the Z^ in- 
variants along the three lines (k yj k z ) = (ir, 7r), (k z ,k x ) = 
(7r,7r), (k x ,k y ) — (7r,7r). We collect them into a vector 
v = {v x ,v y ,v z ) shown in Fig. []](a) (with v x ,y,z — 0, 1). 
It can be shown that v and n together determine the Z^ in- 
variant along all other time-reversal invariant lines. A TSC 
with v 7^ and n = is topologically equivalent to de- 
coupled ID topological superconductor wires aligned in the 
direction of v, each of which has Majorana zero modes at 
the end. Consequently, for any surface plane that is not par- 
allel to v there will be Majorana surface states on that sur- 
face. A generic TSC with both n and v non-vanishing can 
be considered as decoupled layers of 2D chiral TSC coexist- 
ing with decoupled wires of ID TSC. For example, we will 
show later that the RKK model of Sr 2 Ru04 (with the ad- 
ditional spin degeneracy ignored) has the weak topological 
invariants n = (0,0,1), v = (1,1,0), which is equivalent 
to decoupled layers and wires as is illustrated in Fig. [T](b). 
This conceptual decomposition into stacks of lower dimen- 
sional systems will be helpful to illustrate our later discussion 
of dislocation properties. Primary weak topological indices 
were first discussed in the context of time-reversal invariant 
topological insulators |[T0l ITT1 [T8 1 and subsequently both pri- 
mary and secondary indices (and beyond) can be straightfor- 
wardly extracted for the entire periodic table of topological in- 
sulators and superconductors from the K-theory calculation in 
Ref . [13] We briefly remark that although we have defined the 
secondary-index v as a vector the natural structure is actually 
an anti- symmetric two-index tensor which can be interpreted 
as a vector only in 3D. 1191120) 



FIG. 2. (Left panel) Fermi-surface structure of Sr 2 Ru04 showing 
the Fermi surfaces coming from the d xz d yz and d xy bands with or- 
bital mixing. (Right panel) Low-energy quasi-particle spectrum of 
superconducting S^RuCU for a geometry with open boundary con- 
ditions in the ^-direction and periodic boundary conditions on the re- 
direction. The states in red-dashed lines are edge states arising from 
the quasi- ID bands; notice they exist at k x = and k x = it. The 
green dash-dotted lines are the edge states arising from the quasi-2D 
band; notice they exist only around k x = 0. 

Application in Sr^RuO^ We will now discuss the elec- 
tronic structure of the normal metal state of Sr2Ru04 fol- 
lowed by the superconducting pairing scheme as given in the 
paper by RKK[9|. The three relevant orbitals for the elec- 
tronic structure are the t^g multiplet d xz , d yz , d xy which will 
be labelled by a = 1, 2, 3. The layered structure of Sr 2 Ru04 
makes the system behave quasi-two dimensionally; conse- 
quently the first two orbitals are effectively quasi- ID in na- 
ture while d xy is quasi-2D. The bandstructure can be modeled 
using these three orbitals, plus spin, on a simple tetragonal 
lattice with nearest neighbor, and next-nearest neighbor hop- 
pings. The Bloch Hamiltonian, without spin-orbit coupling, 
is 

( e xz {k) A(k) \ 
H(k) = A(k) e yz (k) ®I spin 
\ e xy (k)J 

e xz (k) = —2t cos k x — 2t ± cos k y — 2t\ cos k z 
e yz (k) = —2t cos k y — 2£~ L cos k x — 2t\ cos k z 
e xy (k) = — 2t / (cos k x + cos k y ) — At" cos k x cos k y — 2t\ cos 
A(k) = — 2Asin k x s'mk y 

where values of in-plane hopping parameters taken from RKK 
are t = 1.0, t' = 0.8, ^ = 0.1, t" = 0.3. We have also con- 
sidered an orbital-hybridization term A(k) which arises from 
next-nearest neighbor hopping between different quasi- ID or- 
bitals in xy plane and which removes the crossings of Fermi 
surfaces of quasi- ID orbitals. Hopping amplitudes along z 
are t\ and t\ for quasi- ID and quasi-2D orbitals, respectively. 
Due to the layered structure of the lattice, out-of-plane hop- 
pings are negligibly small |2T] |22) and we shall consider the 
2D limit hereafter. In Fig. [2] we show the Fermi- surfaces due 
to these orbitals. There are three Fermi- surfaces: two around 
(k XJ k y ) = (0,0) and one around (k x ,k y ) = (7r,7r). The 
two quasi- ID Fermi-surfaces from d xz and d yz orbitals do not 
touch as long as A ^ 0. The inner quasi- ID Fermi-surface is a 
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hole pocket, and the outer Fermi- surface is an electron-pocket. 
The round quasi-2D Fermi-surface arises from the d xy orbital 
which we assume is completely decoupled in the 2D limit l23l . 

We now want to consider the properties of the supercon- 
ducting state of Sr2RuC>4. For the quasi-2D band we assume 
triplet p x -\-ip y pairing. For values of the Fermi-level which lie 
within the quasi-2D band (which is expected in experiments) 
this would mean that the system will be a weak 3D topologi- 
cal superconductor with primary index n = (0, 0, 1). For the 
quasi- ID bands we assume the RKK pairing which we now 
describe. Since the quasi- ID and quasi-2D orbitals are as- 
sumed to be approximately decoupled we can write a reduced 
two-orbital model for the quasi- ID orbitals: 

with A = O.lt. The superconducting pairing that RKK pro- 
pose is spin-triplet and intra-orbital. The pairing functions of 
orbital a for this chiral superconducting state are 

A a = id a (k).^, a = 1,2 (3) 
di = zAo sin k x cos k y , (4) 
d2 = iz Ao sin k y cos k x , (5) 

where the direction of d a and the relative phase between 
di and d2 are determined by the small spin-orbit coupling. 
The superconducting pairing winds around the two Fermi sur- 
faces with the same chirality, but since they have the oppo- 
site charge character they contribute oppositely to the wind- 
ing number yielding a vanishing Chern number. However, in 
a clean system there should be edge states located near k x = 
and k x = it if, for example we put the system on a cylinder 
with open boundary conditions in the y-direction and periodic 
boundary conditions in the x-direction. The energy spectrum 
for such a system is shown in Fig. [2] with clear low-energy 
modes near k x = 0, 7r which develop zero modes exactly at 
these k-points. These edge states exist because of a non-trivial 
secondary weak index v = (1,1,0) due to the RKK pairing. 
Even though n, and other quantities, should be doubled when 
the spin degeneracy is taken into account, it has no qualitative 
effect on most of the properties of lattice dislocations we dis- 
cuss below when the effects of spin-orbit coupling are weak. 
We defer the discussions of strong spin-orbit coupling to fu- 
ture work. 

Properties of dislocations and linked dislocation lines: 

In addition to the surface state properties, weak topological 
indices have profound consequences for the properties of crys- 
tal defects such as dislocations! 19, 20ll24l . For the 3D crystal 
we are considering, a lattice dislocation is a line defect around 
which the ions are displaced by an integer- valued vector b in 
the lattice vector basis, which is known as the Burgers vector. 
A dislocation is described by b and another integer-valued 
vector, the tangent vector of the dislocation 1. The relative 
orientation of b and 1 determines the type of dislocation: edge 
(1 • b = 0), screw (1 parallel to b), and mixed (1 neither parallel 
nor perpendicular to b). While b is a topological property of 
a dislocation line, 1 (namely the dislocation-type) is not. 




FIG. 3. (a) Illustration of an edge dislocation line with Burgers vector 
b in a TSC consisting of decoupled layers stacked along the direction 
of topological index vector n. (b) Illustration of an edge dislocation 
line with Burgest vector b and direction 1 in a TSC consisting of de- 
coupled ID wires along the direction of secondary weak topological 
index vector v. The red dots at the end of wires represent Majorana 
zero modes. (c) Illustration of the Majorana zero modes induced by 
linking of two edge dislocations with bi = z and b2 = — x in the 
decoupled plane limit.There is a Majorana zero mode at each dislo- 
cation line, indicated by the red dot and red line. 

Both the primary and secondary weak topological indices 
can be probed by dislocations. The primary weak indices n 
lead to N\ = n • b number of chiral Majorana modes propa- 
gating along the dislocation! 19 ]. Note that N\ is independent 
of the dislocation direction 1 and thus topological. The sign of 
the integer N\ determines the chirality of the edge state, which 
is denned with respect to 1. This fact can be easily understood 
for an edge dislocation with n perpendicular to 1, in which 
case the dislocation can be obtained by adding an additional 
layer of chiral 2D TSC to one side of the dislocation line, as is 
illustrated in Fig. [3] (a). The secondary weak indices v lead to 
non-chiral ID Majorana modes ifiV2 = (bxl)-i/ = l mod 
2. The modes determined by N2 are like the "weak" analog 
of the "strong" modes determined by N\ . The dependence of 
N2 on the variable direction 1 indicates that topological stabil- 
ity will require an additional symmetry which in this case is 
translation symmetry along the dislocation (i.e. the direction 1 
cannot change along the dislocation line). The non-chiral Ma- 
jorana propagating modes are protected by translation sym- 
metry along the dislocation, since its left and right moving 
branches are around the ID momenta and tt, which can- 
not be coupled without breaking translation symmetry. Also 
we see that N2 is nontrivial only for edge dislocations with 
b x 1 / 0. In the topologically equivalent decoupled chain 
limit (which is appropriate for a system with v ^ 0) the dis- 
location bound states can be understood intuitively, as is il- 
lustrated in Fig. [3](b). Decoupled ID Majorana chains along 
the v direction terminate at the dislocation line and the Ma- 
jorana zero modes at their end points couple to form the ID 
non-chiral Majorana edge state. It is thus intuitive to take N\ 
to be akin to a "strong" dislocation invariant and N2 to be a 
"weak" dislocation invariant. 
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One additional effect that has thus far gone unnoticed is 
the property of closed dislocation lines which are linked, i.e., 
dislocation rings. Along a finite-length dislocation ring with 
nontrivial topological invariants N\ and/or N 2 , the Majorana 
fermion energy spectrum becomes discrete, and the boundary 
condition around the ring determines whether there is an exact 
Majorana zero mode at zero energy. Interestingly, the bound- 
ary condition around a dislocation ring depends on its linking 
with other dislocation rings and with flux/vortex lines. To il- 
lustrate this effect, consider the RKK model with v = (1, 1, 0) 
and consider two edge dislocation lines. The first one has 
bi = (1,0,0) and is straight along the z direction so that 
1 = z. The second one is a circle in the xy plane with 
b2 = (0, 0, 1). If these two dislocations are not linked, the 
Majorana fermion mode along the in-plane dislocation has a 
finite size gap with no exact zero mode. This can be under- 
stood easily in the decoupled-layer limit shown in Fig. [3jc), 
in which case the in-plane dislocation circle is the boundary 
of a single-layer disk, and a finite gap of order 1/R (with R 
the radius of the circle) is present. In contrast, when the cir- 
cle encloses the other dislocation line along z direction, in the 
decoupled layer limit the effect is to introduce an edge dislo- 
cation with Burgers vector b x in the disk, which introduces an 
additional Majorana zero mode at the dislocation line. Since 
Majorana zero modes have to come in pairs, there must also 
be a Majorana zero mode around the other dislocation with 
circular shape|25|. Since any generic superconductors can be 
deformed to the decoupled layer limit (due to the absence of 
strong invariant in 3D), the number of Majorana zero modes 
on linking dislocations can be determined only from investi- 
gating the decoupled layer limit, which gives the number of 
Majorana zero modes to be No = Nl (bi x b2) • v mod 2, 
with Nl the linking number of the two edge dislocations. We 
note that the dependence on 1 has dropped out which means 
that TVo is topological and does not require the addition of 
translation symmetry along the dislocations. Thus we see that 
a primary consequence of the secondary weak invariant is the 
determination of bound states on linked dislocations. If we 
had left v as an anti- symmetric tensor this invariant would 
simply be the contraction of the tensor with the Burgers' vec- 
tors of both dislocation lines. 

While the primary weak invariant N\ has no effect in the 
linking of two dislocations, it does determine the Majorana 
zero modes when linking occurs between dislocation lines and 
flux/vortex lines. Obviously when a superconducting vortex 
ring is linked with a dislocation ring, the boundary condition 
for the Majorana fermion along the dislocation will change. 
For odd N\ such a boundary condition change results in a sin- 
gle Majorana zero mode on the dislocation line, and another 
one on the vortex line. The exitance of Majorana mode on the 
dislocation line is determined by Nq = NlNi mod 2, where 
Nl is the linking number between a dislocation line and a 
vortex line. 

Physical consequences: Before we conclude, we now dis- 
cuss a few measurable experimental consequences of the pre- 
dicted Majorana fermions bound to dislocations: 
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FIG. 4. Spatial profile of the lowest energy bound state that is bound 
to the ends of a dislocation that stretches from x = 40 to x = 120 
on the line y = 50. The total lattice size is 160 x 100. The asym- 
metry in the density profile is due to finite- size error introduced by 
discretization. 

1) The one-dimensional Majorana modes have a constant 
density of states. Consequently, it contributes a specific heat 
Cy that is linear in temperature, and is proportional to the den- 
sity of dislocations 1 26 ]. Since the dislocations in a Sr2Ru04 
sample can be observed by transmission electron microscopy 
(TEM)[ 27 ], measuring samples with different dislocation den- 
sity can verify whether Cy /T in the low temperature limit is 
indeed proportional to the dislocation density. 

2) Due to the nontrivial primary weak topological invari- 
ant n, the dislocations with Burgers vector b = z have chiral 
Majorana fermion modes. Thus, each such dislocation carries 
a chiral heat current Ie = ^I/T lf28l along the dislocation 
line. With a random distribution of dislocations in the system, 
the flow direction of the chiral heat current is also random, 
so that a net chiral heat flow will not be observed. How- 
ever, the chiral heat current along random dislocations will 
contribute a thermal conductivity that is proportional to the 
dislocation density. Furthermore, it is possible to have a sys- 
tem with imbalanced dislocation lines. For example, in a sys- 
tem with a uniaxial strain u zz (y) that has a gradient along the 
^-direction, it is preferred to have edge dislocations along the 
x-direction with more Burgers vector b = +z than b = — z. 
If the density of the b = +z (b = — z) edge dislocations are 
p + (p_) respectively, there will be a chiral heat current den- 
sity ]e = (p+ — p~) ^24^ • Compared to the thermal cur- 
rent carried by the edge states which is easily overwhelmed 
by bulk thermal conduction, the dislocation current can be a 
bulk effect that remains finite in the thermodynamic limit. 

3) If there is an edge dislocation ring in the xy -plane (with 
Burgers vector bi = z) and a second edge dislocation along 
the z-axis threading the ring (with Burgers vector h 2 in- 
plane), then the Majorana zero mode in the second dislocation 
may be observable by scanning tunneling microscopy (STM). 
In particular, when the first dislocation is at a crystal surface, 
it will be a disk-shaped plateau on the surface, threaded by 
the second dislocation line and be easy to locate using, for 
example, TEM. 
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